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Abstract 



We study the statistics of level widths of a quantum dot with extended con- 
tacts in the absence of time-reversal symmetry. The widths are determined by 
the amplitude of the wavefunction averaged over the contact area. The dis- 
tribution function of level widths for a two-point contact is evaluated exactly. 
The distribution resembles closely the result obtained when the wavefunction 
fluctuates independently at each point, but differs from the one-point case. 
Analytical calculations and numerical simulations show that the distribution 
for many-point contacts has a power-law behavior at small level widths. The 
exponent is given by the number of points in the lead and diverges in the 
continuous limit. The distribution of level widths is used to determine the 
distribution of conductance peaks in the resonance regime. At intermediate 
temperatures, we find that the distribution tends to normal and fluctuations 
in the height of the peaks are suppressed as the lead size is increased. 
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I. INTRODUCTION 



Usually any measurement of conductance in a metallic system assumes at least two 
attached leads. The existence of these leads connecting the system to reservoirs broadens 
the electronic energy levels. For a macroscopic sample, the width of an energy level is 
typically much larger than the distance between neighboring levels. As a result, one observes 
a smooth dependence of the conductance on the Fermi energy. 

Recently, the advances in nanometer technology have made possible the fabrication of 
very small semiconductor devices, known generically as quantum dots,0^ where one has a 
fixed number of conduction electrons confined into an island. In these systems, one can 
narrow the level width by reducing either the size of the leads or their transmittance. In 
the first case it is customary to speak about channels: If we consider an ideal lead as a 
wave guide with a given cross section, we can associate a channel to each state due to 
transverse quantization. The number of channels is approximately the area, in units of 
electron wavelength, of the contact between lead and dot and each channel is characterized 
by its transmittance. 

The conductance of a system can be calculated through the well-known Landauer- 
Biittiker formula.lii Following this approach, the distribution of conductances of a quantum 
dot has been evaluated for two distinct cases: for weakly coupled, pointlike, leadsi of a 
closed system, and for leads with any number of channels attached to a ballistic cavity.li 

Here we will consider the situation at low transmittance, when leads behave as tunnel 
contacts. In this case, each electron eigenstate corresponds to a peak in the Fermi energy 
(gate voltage) dependence of the conductance.Bi This allows one to make a real spectroscopy 
of electrons in quantum dots. The height of each peak is determined by the probabilities 
of tunneling through the leads, as well as by the amplitude of the wavefunction near the 
contacts. The later means that these heights are randomly distributed. The distribution 
function of the peaks has been determined for the case of pointlike contacts in Refs. and 



T0| . In this paper we will deduce the distribution considering leads of arbitrary size. 

We begin by studying the statistics of level widths, which is connected to the fluctu- 
ations of conductance peaks through the Landauer-Biittiker formula. The supersymmetry 
technique and the nonlinear a mo delll3§ provide the framework for the exact calculation 
of the distribution of level widths for leads with two-point contact. The distribution dif- 
fers drastically from the single-point case, but we find that inclusion of correlation between 
the wavefunction fluctuations at the two points in the lead does not significantly alter the 
overall form of the distribution of level widths. This analytical result is confirmed through 
numerical simulations of a quantum dot. The exact distribution of conductance peaks is also 
evaluated for two-point leads and we find that it deviates from the single-point distribution 
mainly for small peak heights, where it shows a linear dependence. 

For leads with many-point contacts we do not know of any method which enables a 
general evaluation of the distribution function of level widths. Thus we assume that the 
wavefunction at each point fluctuates independently and proceed to calculate the total dis- 
tribution. By comparing numerical with analytical results, we show that this approximation 
yields a good qualitative understanding of the large- lead limit, since correlations do not 
effect the behavior of the distribution very strongly. 

An important use of quantum dots is found in the study of chaos.lii"Ei Because one can 
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obtain extremely clean samples and shape them into different forms, it is now possible to 
fabricate the so-called quantum billiards, where one has a unique chance to study quantum 
and semiclassical physics. Although the supersymmetry method assumes averaging over 
disorder, one can conjecture on very firm grounds (the ergodic hypothesis) that it applies 
quite generally to quantum chaos problems. Therefore, we expect that our results should 
be able to describe any system where the underlying dynamics is chaotic, regardless as to 
whether it is diffusive or ballistic 

The organization of this paper goes as follows. In Sec. we describe the model used and 
the basic concepts related to conductance peaks and level widths. The way the numerical 
simulations were done is explained in Sec. |T|. Some of the results found in Ref. ^, which 
initially motivated the present work, are reviewed in Sec. The case of two-point leads 
is treated exactly in Sec. [V| and the case of many-point leads, under the approximation 
of independent point fluctuations, is left to Sec. 0. Finally, in Sec. |VI1| we draw our 
conclusions and point out some experimental consequences of our work. 



II. FORMULATION OF THE PROBLEM 

The electronic system we consider is formed by noninteracting spinless fermions confined 
to a small region through some strong confining potential. The system is probed by two 
leads, described by discrete sets of points, which are weakly coupled to the system. More 
precisely, we have the Hamiltoniani' ^ 



H = -— + U{r) + 



2m ' ' 2'KAf 



(1) 



where U{r) is the potential (containing the effect of random impurities and confining walls) 
and M = 1/ (V^A) is the mean density of states, with V denoting the total volume and A the 
average level spacing. The term inside the brackets represents the contact leads, denoted by 
2, and ai^2 are the dimensionless coupling parameters. 

In a situation characteristic of scattering experiments, when ai^2 ^ 1 and the leads are 
open, the number of channels is equal to the number of wavelengths which fit into the lead 
cross section. As we shall demonstrate, for the weakly coupled leads of a closed quantum 
dot, the distinction between pointlike (single channel) and multichannel leads is more subtle 
and deserves a careful analysis. 

When tunneling through the leads in totally suppressed (ai,2 = 0) the energy levels have 
zero intrinsic width. The levels will gain a small but finite width if the coupling to the 
leads is made very weak, i.e., if < ai^2 ^ 1. As one can show by first order perturbation 
theory, the broadening will be caused by local fluctuations of the electron wavefunction in 
the regions of contact with the leads. Calling the total level width Ty, we have 

= 7i/,i + lu,2 , (2) 

with the partial level widths given by 



^)EIV^^WI' (^ = 1,2), (3) 
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where ipi, is a single-particle eigenfunction of the Hamiltonian in Eq. (|l|) for ai 2 = 0. 

The knowledge of the statistical properties of 7^ is crucial to describe some important and 
measurable effects in quantum dots. From the distribution of level widths we can extract the 
distribution of conductance peaks in the resonance regime. It is rather simple to understand 
why: Let us assume that the conductance at a given energy E can be evaluated through the 
Landauer-Biittiker formula,!'! 



(4) 



where e^, is the eigenvalue associated with the eigenfunction ipy. In the resonance regime 
(7, T <^ A), only the eigenstate whose energy is the closest to E significantly contributes to 
the conductance. G{E) is very small when E is between adjacent energy levels, and grows 
rapidly when E — > Eu- Therefore, at a peak of the conductance, we have 

^^ = ¥7^^ ^ ^ \Mr.)\'\Mr.)r (5) 
Using the last equation, plus Eqs. (|D and (H), we obtain the Breit-Wigner formula 
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which is correct at zero temperature, or when T -C 01,2^. In order to describe a wider range 
of temperatures, we have to use instead the well-known relation 

2e^ r de lu,ilu,2 , -2 f ^ 



Let be the number of points in each lead. When the leads are placed very far apart, 
the connection between the distribution of level widths, PN{lu,i), and the distribution of 
conductance peaks, RN{gm), is given by the convolution 

RN{gm) = d-iu^PNinv,!) c?7i^,2^7v(7i/,2) ^ i^Q^ - G^h/2e^j . (8) 

From Eq. (|^) we can extract another limit, which occurs when the temperature exceeds 
the intrinsic energy level width, T ^ 0:1,2 A, but is still low enough so that there are well- 
resolved resonance peaks (T <C A). In fact, this intermediate regime is typical of some 
experimentslll and is described by the Hauser-Feshbach formula 

" h l2Tj7..i + 7.,2 ■ ^ ^ 

In practice, one can generate many peaks in the conductance by varying the energy E 
(i.e., by sweeping the gate voltage in the quantum dot), or, else, by varying some external 
parameter, say, an applied magnetic field. Similarly, in numerical calculations one can apply 
a Aharonov-Bohm flux to the system, and then generate a family of curves e^{(f)), which 
yields a random sequence of conductance peaks, G^{(j)). As long as the system is ergodic, all 
these procedures are equivalent to the construction of an ensemble of different realizations 
of disorder (in the case of a diffusive regime) or boundaries (if the regime is ballistic). 
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III. NUMERICAL SIMULATION 



Before we move to the analytical calculations, let us describe the numerical simulation 
of a disordered quantum dot which we used to illustrate the main points of our work. The 
numerical results will be shown along the discussion of the analytical ones. Our motivation 
was to obtain the histogram PnIi) as a function of the separation between the points within 
the lead. We used a two-dimensional Anderson model of noninteracting spinless electrons 
with nearest-neighbor hopping and diagonal disorder, i.e., the Hamiltonian 

N 

Ha = - ^{c\cj + c]ci) + ^Wi c\ci , (10) 

with Wi uniformly distributed in the interval [—W/2 ,W/2]. We restricted our calculations 
to only one realization of {wi}; hence, no average over disorder was performed. On the 
other hand, we imposed quasiperiodic boundary conditions to the electrons by making the 
geometry toroidal and introducing Aharonov-Bohm phases in the hops, 

clc, cfc.e-'^'^ , (11) 

where (pij assumes one of two values, (px or (py, depending on whether the hop occurs along 
the X or y axis, respectively. By varying (p^ and (py within the interval (0, vr) we obtained a 
large set of states, thus creating the statistical ensemble needed to construct the histograms. 

For a given set of phases (px and (py the Hamiltonian was diagonalized through standard 
procedures and all eigenstates obtained. The leads were then chosen as sets of sites placed 
throughout the grid, such that we could probe the density fluctuations at different regions 
of the dot. Notice that the leads were passive objects, which did not affect the eigenstates. 

Briefly, we mention the existence of a technical complication, intrinsic to simulations of 
quasiperiodic systems, which is the crossover with respect to time-reversal breaking. The 
threshold for breaking time-reversal symmetry with Aharonov-Bohm fluxes is proportional 
to the conductance of the system.0 If the disorder is too strong or the system too small, a 
complete time-reversal symmetry breaking may never be achieved due to the periodic nature 
of fluxes. As a result, in order to fall into the class of the unitary ensemble and be able 
to compare the simulations with the calculations to follow, we needed to work with a large 
grid (32x23) and keep the disorder very low. The system was then maintained marginally 
diffusive (W ^ 1). 



IV. POINTLIKE LEADS 

We will consider first the simple case of pointlike leads, namely, when the typical size of 
the contact is smaller than the area A'^"^, where A is the electron wavelength. In the lattice 
version of the problem, this certainly happens when the lead consists of a single site. For a 
chaotic system with broken time-reversal symmetry, the amplitude v = V\il!y{r)\'^ fiuctuates 
according to an exponential distribution. This result follows from the assumption that the 
components of the eigenfunction are not correlated among themselves and are uniformly 
distributed.ii Therefore, by connecting v to the level width 7^ [see Eq. (^] we can write 
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(To simplify the notation, we will hereafter drop the eigenstate label.) In Fig. 1 we have 
plotted Pii'ji) obtained from our numerical simulation against Eq. (|12D for two different 
regions of the spectrum. The states around the bottom of the band used in the evaluation 
of the level width distribution were checked to be extended over the entire grid. 

If the distance between the leads Ai and A2 is much larger than the electron wavelength 
A (so that the wavefunction fluctuations at ri and r2 are independent), we can use Eq. 
(I) or Eq. (H), together with Eqs. (|]) and (0), to evaluate analytically the distribution of 
conductance peaks for single-point leads, -Ri (gm) ■ This calculation has been recently done in 
the literature.ii In Ref. |] the authors found an expression for the conductance distribution 
at intermediate temperatures not only for single- channel leads, but for multichannel ones as 
well. Their derivation, based on random matrix theory, enabled them to consider the case of 
time-reversal symmetry, although they were constrained to treat only the case of independent 
channels. In addition, for simplicity, they treated only the situation of symmetric leads. A 
similar work is found in Ref. where the authors have also studied the effect of chaotic 
and regular dynamics in the distribution of conductance peaks. 

On the other hand, the authors in Ref. ^ performed their analysis exclusively for the 
case of broken time-reversal symmetry and pointlike leads, but took into account a possible 
asymmetry between the leads. Since we shall extend their work, we will display and comment 
on some of their results. For T ^ 0:1,2^ they obtained the following distribution 

P , . _ ejl-gj [l + (2-^^)(a^-l)] 
where the asymmetry parameter a is given by 

. (14) 

Notice that, according to Eq. (p!^ , for a = 1 (symmetric leads), the most probable value 
of gm is 1. However, for a > 1 the distribution begins shifting to lower values of gm- For 
a ^ 1, the most probable value of gm tends to zero. This is exactly what one would expect 
to happen in a resonance tunneling experiment at zero temperature: The conductance is 
maximum when the barriers are identical. 

For the intermediate regime, when ai^2^ ^ T <^ A, it is convenient to rescale the 
distribution to obtain! 

^^^^^ ^ ^/^^R^^gj = a:e-^^[K,{x) + aK,{x)] , (15) 
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where 



X = ^-^^"^ (161 

and Kn{x) is the modified Bessel function of order n. Notice that the rescaling of gm to 
X and of Ri to TZi turned the distribution into a temperature-independent function, whose 
form depends only on the asymmetry between the leads. We also point out that for both 
temperature regimes, one has an finite offset in the distribution, i.e., Ri{gm — > 0) 7^ 0. 
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V. TWO-POINT LEADS 



In the previous section we reviewed known results. In this and the following sections the 
main goal will be to study the effect of extended leads in the distribution of level widths 
and conductance peaks. We start by analyzing the situation of two-point leads, where the 
level width is given by 

l^='^{v^+V2) . (17) 
vr 

The amplitudes are defined as Vi = V\ilj^{ri)\'^, i = 1,2. When the points 1 and 2 are far 
apart (r = |ri — ^ A), the amplitudes fluctuate independently and the distribution of 
level widths is equal to the convolution of the distributions for the isolated amphtudes [see 
Eq. ([121)]. Consequently, 

P.(7.)^(^)%.exp(-||). (18) 

Notice that the most probable value of ji is not zero, as for the one-point case, but ajA/vr. 
This fact by itself signals a strong qualitative change in the distribution. What happens, 
then, when we move the points closer together, so that r ^ A? This question can only be 
fully answered if we know the joint probability distribution 

Q2{vuV2;r) = {6{v^-V\Mri)\') S{v2 - V\Mr2)\')) ■ (19) 

However, a simple analysis can be made before one starts seeking the exact form of 
Q2ivi,V2;r): Since there is no special reason for Q2(0,0;r) to be singular if r 7^ 0, we 
should have ^2(71 — >■ 0) oc 7,. On the other hand, (^2(0, 0; 0) is apparently singular in order 
to Eq. (P) hold true. 



The exact calculation of Q2(^i, ^2; 1^) can be performed by the supersymmetry method.E2l 
The result depends on the symmetry of the Hamiltonian; here we will consider only the case 
of broken time-reversal invariance (unitary ensemble), which yields 

Q^ivuV2;r) = ^-1^ exp [-^^^) [j^^^V^] ' (20) 
where /(r) is the Friedel-like function 



Jo(27rr/A) {d = 2) 

(21) 

(A/27rr) sin(27rr/A) [d = 3) , 

and Io{x) and Jo{x) are Bessel functions of order zero. The exact expression for the two-point 
distribution of level widths can be readily calculated: 
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-P2(7i) = / dvi f dv2 Q{vi,V2;r) 6{-fi - {aiA/7T){vi + t;2)) 



TT 



a^Af{r) 




The limit r ^ A (points placed far apart) can be easily extracted from Eq. (|22|) and it 



agrees with Eq. (|T^. For the opposite limit, r = 0, when the points coincide, one gets 

^■.(7.) = {^] exp {-^) . (23) 



which is the distribution for a pointlike lead with a coupling constant twice as large, as we 
would expect. 

The large fluctuation tail of the distribution tends to an exponential function, just as 
Eq. (0). The correlation between points, however, renormalizes the coupling constant Ofj 
to + /(r)]. Notice also the linear behavior at 7^ — > 0, for any r 7^ 0, consistent with 
what we anticipated. The slope of the distribution at small widths becomes steeper as the 
distance between points decreases, but the qualitative aspect is rather independent of /(r) 
and one can usually approximate P2{,li) by Eq. (^). 

We stress that the main effect of point correlation within the lead is that the characteristic 
number of channels of a lead is not quite accurately given by the cross section divided by the 
electron wavelength. We can say that there is a crossover between single and double-channel 
behavior driven by the distance (correlation) between the points, i.e., by the parameter /(r). 

In Fig. 2 we show some distributions of level widths for two-site leads with different 
site separations obtained from our numerical simulation. One observes that the correlation 
between sites is small, even for r = 1 (in lattice units) and at the low energy portion of the 
spectrum (where A is larger than at the middle of the spectrum). This effect, peculiar to the 
tight-binding model adopted here, can also be visualized in Fig. 3, where we have plotted 
the density autocorrelator as a function of site separation. 

Given Eq. ( P^ we can proceed to evaluate R2{gm)- In analogy to the pointlike case, 
an analytical treatment is viable only for the two temperature regimes T -C q;i,2^ and 
ai,2A <C T <C A; for T comparable to Q!i,2A, one has to compute the distribution numerically 
and the dependence on the temperature cannot be rescaled out of the distribution. For 
simplicity, we will concentrate on the analytical results. We first present the expression we 
obtained for the resonance conductance distribution in the T <^ ai 2 A regime: 



R2{gr. 



9{l-gJ ([l + /(^)2][i + (2_^J(a2_l)] 

(24) 



Af{rf^T^ \ [i+^^(a2-i)]2 

[1 - /(r)2][l + (2 - g„,){al - 1)] [1 - /(r)2][l + (2 - g„.){al - 1)] 



2[l + gm{a\ - l)f 2[l + g^{al - I)] 



where 



(a±/(r)v/^^)/^l-/(r)2, (25) 

and /(r) is defined in Eq. (^1]). In Fig. 4 we show R2{gm) for different values of /(r) and 
the asymmetry parameter a. Notice that the most probable value of gm is maximum when 
the leads are symmetric, as we have argued in the previous section. 



8 



Clearly, by setting r = in Eq. p^ ) we recover Eq. (0). In the opposite limit, r ^ A, 
we can expand all terms to 0{p) and obtain 

Mdm) = -—==-— — -— . 26 

A y/1 - gm [I + gm[a^ " 1)J* 

The distribution of conductance peaks in the regime ai^^^ ^ T ^ A can again be 
expressed as a universal, temperature independent function if one uses the rescaling shown 
in Eqs. (|T5|) and (|16|) . After some straightforward manipulation we obtain 

^2(^) = j:{-ire-^^''^[K,{x,) + a,K,{x,)] , (27) 

^Iv) p=o 

where Xq = x/[l + /(r)], X2 = x/[l - /(r)], Xi,3 = x/y'l - /(r)2, ao,2 = a, Oi = a+, and 
a3 = a_. In Fig. 5 we have plotted 'R-2{x) for different values of /(r) and a. When we set 
r = in Eq. (|27|) , we recover Eq. (^) with coupling constants twice as large. On the other 
hand, in the limit r ^ A, one has, after an expansion to 0(/^), 



7^2(a;) = — — 



x{a^ + l)iro(a;) + (2ax + 20^ - . (28) 



In both temperature regimes, the most notable difference between R2{gm) and Ri{gm) is 
in the small gm behavior: R2{gm 0) is linear in gm for any r 7^ 0, in contrast to Ri{gm), 
which tends to a constant in the same limit. Correlations between points within the leads 
do not affect the linear dependence. 



VI. iV-POINT LEADS 

Some of the conclusions we have drawn for the two-point lead are straightforward to 
extend to the A^-point case. For instance, let us call Qn{vi,V2, ■■■,viy) the joint probability 
distribution of the amplitudes Vk = V\ip{rk)\'^ at A^ points. Then, since Q^r(0, 0, 0) is 
nonsingular for any arrangement where points do not coincide (i.e., when no two points are 
completely correlated), we always have Pni'^i — ^) oc 7^^"^ 

Unfortunately, the supersymmetry technique becomes impractical to use in any deriva- 
tion where N > 2. We do not know of any other method suitable for this task either. 
Therefore, in this section we will simply assume that the point fluctuations within the leads 
are completely independent, which is certainly correct for small enough values of 7^ or gm- 
We will then derive expressions for PNili) and RN{gm) and look at the numerical simulations 
to gain some insight about the general case, when correlations can be present. 

If the amplitudes Vk fluctuate independently, we are allowed to convolute A^ distributions 
given by Eq. (|1^) to obtain the distribution of level widths for a A^-point lead. In this way, 
we find that 
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where 7j = (ajA/vr) X^^Li "^^/o and, according to our definition, (7j) = NaiA/n. In the hmit 
iV ^ 1 one can check that the distribution becomes Gaussian around 7^ ~ OiAN/iT, with a 
width proportional to 1/ VN, 



The results of our simulations for =4, 9 and 16 are shown in Fig. 6 for two different 
lead geometries and for extended eigenstates around the bottom of the energy band. The 
deviation from Eq. (|30| ) becomes substantial for large N when the lead sites are close 
together. This indicates that the modest site-to-site correlation (see Fig. 3) can add up 
to make the effective number of channels smaller than the number of sites, particularly for 
small widths. As we move the sites farther apart, correlations become negligible and the 
distributions start to agree with Eq. (|30D . 

From Eq. (p9|), together with Eqs. (P), (^, and (||), we can proceed to evaluate RN{gm)- 
In the very low temperature regime, T -C 01,2^, we obtain 

N{9m) 2r(^)^/r^[i + ^^(^2_i)]2iv A. (2N-2k)\{2N + 2k)\ ' ^ ' 

which contains all the essential features of Ri and R2- In the other regime, where ai_2A ^ 
T -C A, we obtain, after an appropriate rescaling. 



i<^n{x) - ^-j r(iv)2 [jmr h {N-i)\{N + i)\ / • 

One can easily check that Eq. (^) becomes equivalent to Eqs. (|T^) and (|26|) for = 1 and 
N = 2, respectively. Similarly, Eq. (^) is reduced to Eq. ( [T3| ) for A^ = 1, and to Eq. ( PBD for 
N = 2. We remark that both Eqs. (|3TD and (^) yield the same limit RN{gm 0) oc g^'^y 
i.e., the distributions behave like a power law for small peak heights. 

The large A^ limit makes the distribution of conductance peaks narrower. For instance, 
when A^ 1, Eq. ( p2D tends to a Gaussian law, centered around x ~ 0{N), 



2^-1 (2Ar)!e-'^- j Ko{x) , [(«i/«2y/' + («2/«i)'/']i^K^) 



a ( ix-2N/af\ , , 



Because {gm) ~ 0{N), but 6gm ~ 0{yN), fluctuations in the height of resonance peaks 
are suppressed for very large leads. This should be contrasted to the universal fluctuations 
of conductance as a whole, which obey {gm)l5gm ~ Oil). Furthermore, {g^) oc (1/ai + 
1/02)"^ 5 which is the classical result (Ohm's law). 



VII. CONCLUSIONS 

We have considered the statistics of the widths of electronic states inside a quantum 
dot coupled to a bulk reservoir through contact leads. It is well-known that in the case of 
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pointlike contacts the distribution of widths obeys a universal exponential law, which corre- 
sponds to Gaussian fluctuations of the wavefunction. We have shown that the statistics of 
the widths is strongly influenced by the size of the contact area between the leads and the 
dot. For two-point contacts we presented the exact expression for the distribution function 
of widths and found that it closely resembles the result one obtains by assuming no corre- 
lation between the points. The correlation effects appear only for large fluctuations, whose 
contribution to any observable quantity is exponentially small. Following this observation, 
we suggested that the A^-point lead can be considered as a set of N points fluctuating in- 
dependently. We evaluated the distribution of widths analytically and compared the results 
with direct measurements from numerical simulations. A good qualitative agreement was 
found. 

The central region of the distribution tends to be Gaussian as N increases. Therefore, 
the fluctuations in the level width are suppressed for large leads and we can describe the 
dot in terms of the average parameters. In the region of very large widths, the distribution 
does not follow a Gaussian law, but rather decays exponentially. For small widths, the 
distribution behaves as a power law. 

For the continuous system with extended leads we can expect a similar behavior. The 
distribution of widths is Gaussian around the average, which is proportional to the effective 
number of channels A^^o = A/A°'~^, where A is the cross section of the lead contact. The tail 
of the distribution for very large widths is exponential. In the region of small widths the 
distribution behaves as 

P(7.) (X 7f . (34) 

The effective number of channels A^ is equal to the number of independent "points" in the 
lead, which increases for small 7^ as A^ cx A^o/7j- 

The relevance of these results to experiments appears in the determination of the distri- 
bution of conductance peaks in the resonance regime. The main consequence of the finite 
extension of the leads is found in the power law behavior of the distribution at small peak 
heights. As the lead becomes very large, the distribution tends to a Gaussian law. The 
possibility of asymmetry between the leads is taken into account in our work. Finally, we 



point out that Eqs. (24) and (|27|) can be used to describe the statistics of four-probe dots 



where the constrictions are made very narrow (pointlike leads). 
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FIGURES 



FIG. 1. Histogram of level widths for single-site leads at two regions of the energy band: middle 
(circles) and bottom (squares). The dashed line is the universal prediction (exponential law). The 
level width is rescaled so that (71) = 1. 

FIG. 2. Histograms of level widths for two-site leads at two regions of the energy band: middle 
(squares) and bottom (circles). The level width is rescaled so that (72) = 2. r is the distance in 
lattice units between sites within the lead. The solid and dashed curves are the predictions for 
/(r)=0.35 and /(r) = 0, respectively. 

FIG. 3. Density-density correlator as a function of site separation for two regions of the energy 
band: middle (circles) and bottom (squares). The density is defined as p(r) = |V'i^(r)p and the 
average is performed over v (eigenstate) and grid location. 

FIG. 4. Distribution of resonance conductance peaks in a dot with two-point leads at very low 
temperatures (T <C ai,2A). The main plot shows the curves for the symmetric case (a=l) and 
different separation of lead sites. Notice the liner behavior as 5^ — for any fT^l. The insert 
shows the dependence with the asymmetry parameter for a fixed distance between sites (f=0.5). 

FIG. 5. Distribution of resonance conductance peaks in a dot with two-point leads at in- 
termediate temperatures (01^2^ ^ T ^ A). The conductance is expressed in the rescaled form 
X = ^Tgm/ \/a\a2^- The main plot shows the curves at a=l (symmetric case) and variable dis- 
tance between the lead sites. The insert shows the dependence on the asymmetry parameter for a 
fixed f=0.5. 

FIG. 6. Histogram of level widths for A^-site leads and for eigenstates at the bottom of the 
energy band. The level width is rescaled so that (7^) = N . The dashed line is the prediction for 
independent site fluctuations of the wavefunction. The squares correspond to leads where sites 
were placed together (r = 1) and the circles correspond to leads with r = 3. As an example, in the 
insert we show how a iV=4 lead with r = 2 was implemented. 
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